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ABSTRACT 



This paper addresses the origin of the far-infrared (FIR) continuum of QSOs, 
based on the Quasar and ULIRG Evolution Study (QUEST) of nearby QSOs and 
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ULIRGs using observations with the Spitzer Space Telescope. For 27 Palomar- 
Green QSOs at z<0.3, we derive luminosities of diagnostic hues ([Nell] 12.8/im, 
[Nev] 14.3;um, [O iv] 25.9/im) and emission features (PAH7.7/im emission which 
is related to star formation), as well as continuum luminosities over a range of 
mid- to far-infrared wavelengths between 6 and 60;um. We detect star- format ion 
related PAH emission in 11/26 QSOs and fine-structure line emission in all of 
them, often in multiple lines. The detection of PAHs in the average spectrum of 
sources which lack individual PAH detections provides further evidence for the 
widespread presence of PAHs in QSOs. 

Similar PAH/FIR and [NeH]/FIR ratios are found in QSOs and in starburst- 
dominated ULIRGs and lower luminosity starbursts. We conclude that the typ- 
ical QSO in our sample has at least 30% but likely most of the far-infrared 
luminosity (~ 1O^'^'"^^L0) arising from star formation, with a tendency for larger 
star formation contribution at the largest FIR luminosities. 

In the QSO sample, we find correlations between most of the quantities stud- 
ied including combinations of AGN tracers and starburst tracers. The common 
scahng of AGN and starburst luminosities (and fluxes) is evidence for a starburst- 
AGN connection in luminous AGN. Strong correlations of far-infrared continuum 
and starburst related qiiantitics (PAH, low excitation [Nell]) offer additional sup- 
port for the starburst origin of far-infrared emission. 

Subject headings: Infrared: galaxies - Galaxies: active - IR observations - 
Galaxies: starburst - quasars: emission lines 

1. Introduction 

The infrared properties of luminous active galactic nuclei (AGN) and ultraluminous 
infrared galaxies (ULIRGs) hold clues for the understanding of galaxy formation, the star 
formation history of the universe and the connection between black hole and galaxy for- 
mation and evolution. Physical insight into the co-evolution of AGN and star formation 
can be gained by study of the low redshift members of these populations. The Infrared 
Astronomical Satellite (IRAS) and the Infrared Space Observatory (ISO) have allowed the 
flrst comprehensive far-infrared (FIR) studies of extragalactic sources and some high quality 
mid-infrared (MIR) observations on the brighter sources (see e.g. the review of Genzel & 
Ccsarsky 2000). Issues discussed in IRAS and ISO studies include the nature of ultralu- 
minous infrared galaxies and their starburst and AGN energy sources (Genzel et al. 1998; 
Klaas et al. 2001), and the SEDs of QSOs and their evolutionary implications (Haas et al. 
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2003). A 2.5 to 45 //m spectral inventory of starburst and AGN prototypes (Sturm et al. 
2000) and detailed mid-infrared spectroscopic investigations of local lower luminosity AGN 
(Clavel et al. 2000; Sturm et al. 2002, hereafter S02) have addressed AGN unification and 
the starburst- AGN connection. 

The Spitzer Space Telescope (hereafter Spitzer) has enabled us to build upon previous 
work with significantly improved sensitivity. This work is part of a series of papers describing 
the results of the Quasar and Ulirg Evolution STudy (QUEST), focussing on mid-infrared 
spectroscopy of a sample of 54 QSOs and ULlRGs at redshifts z^O.3 with the Infrared 
Spectrograph IRS onboard Spitzer. The sample was selected with the aim of investigating 
the possible connections between those two groups of luminous (LboI ^ lO^^I/©) objects 
in the nearby universe, and is closely connected to optical and near-infrared studies of the 
morphology and dynamical properties of these two populations (e.g. Dasyra et al. 2006a,b; 
Veilleux et al. 2006). The purpose of the present paper is to study starburst signatures 
in QSO host galaxies and to show that most, and in some cases all, of the far-infrared 
luminosity in QSOs showing strong PAH features is due to starburst activity. A forthcoming 
complementary paper addresses the spectral energy distribution (SED) of QSOs including 
the clear far-infrared starburst contribution to the A > 30/xm continuum of QSOs. 

To investigate the link between AGN activity and star formation and the extent to 
which they occur simultaneously, it is important to quantify the star formation activity in 
QSO hosts. Such measurements arc made difficult, however, by the observational problems 
of detecting star formation tracers in the presence of extremely powerful AGN emission. 
SED studies based on the IRAS and ISO space missions have established QSOs as sources of 
(sometimes) strong far- infrared emission, (e.g. Neugebauer et al. 1986; Haas et al. 2003). In 
addition to a nonthermal continuum that is detectable in the infrared only in flat spectrum 
radio- loud QSOs (e.g. Haas et al. 1998), a strong far-infrared emission component is often 
observed, at varying levels with respect to the strong AGN mid-infrared continuum. Due 
to its steep falloff in the submillimeter regime, the origin of this far-infrared emission must 
be thermal emission of optically thin dust (Chini et al. 1989; Hughes et al. 1993). While 
the warmer T~200K dust, which dominates the mid-infrared SEDs of QSOs, is generally 
accepted to be predominantly AGN heated, there is still considerable dispute about the 
origin of the cooler T~50K emission often dominating the far-infrared. Direct heating by 
the powerful AGN, but at distances ensuring sufficently low temperatures, is one possibility 
(e.g. Sanders et al. 1989; Haas et al. 2003; Ho 2005). Other models prefer an origin in vigorous 
star formation in the QSO host (e.g. Rowan-Robinson 1995). Rowan-Robinson (1995) used 
radiative transfer modelling to infer an SED of AGN heated dust that, in vL^, units, peaks in 
the mid-IR and decays towards the far-infrared, a feature shared by many other such models. 
In the QSO SEDs that are often flat over a wide wavelength range including the far-infrared. 
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the far-infrared component is then plausibly ascribed to a component with an SED similar to 

that of a star-forming galaxy, in accordance with evidence for coexistence of star formation 
and AGN in spatially resolved lower luminosity AGN. Rowan- Robinson (1995) found a tight 
correlation of AGN optical emission and mid-infrared continuum and a weaker correlation 
between optical and far-infrared emission, which is supporting the view that the far-infrared 
does not result directly from AGN heating but that there is a connection between AGN and 
starburst luminosities in the QSOs. 

Our goals are to quantify star formation in QSO hosts and to estimate its contribution 
to the the far-infrared emission. In the mid-infrared, the contrast between the emission 
from possibly dust-obscured star formation and from the central AGN is favourable, and 
estabhshed star formation tracers are available. We use two such tracers: 

(1) The mid-infrared broad aromatic 'PAH' emission features arise in regions of the 
interstellar medium of a galaxy where their aromatic carriers are present, and where their 
transient excitation is made possible by a non-ionizing (< 13.6eV) soft UV radiation field. 
This is the case in the photodissociation regions (PDRs) that accompany Galactic star 
formation regions (e.g. Verstraete et al 1996), as well as in the diffuse interstellar medium 
where they are excited by the general interstellar UV radiation field that has leaked from 
its OB star origins to large scales (e.g. Mattila et al. 1996). PAHs have been used as a 
quantitative tracer of star formation in galaxies (e.g. Genzel et al. 1998; Forster Schreiber 
et al. 2004; Calzetti et al. 2005). Metallicity above 0.2 solar is a prerequisite for strong PAH 
emission (Engelbracht et al. 2005), a condition that is probably safely met for local QSO 
hosts. Destruction of the PAH carriers in energetic environments but survival in starburst 
PDRs (though not in H ll regions proper) is key for the use of PAH features as diagnostic. 
The PAH features are absent from the hard radiation environment of AGNs according to 
both empirical (e.g. Roche et al. 1991; Le Floc'h et al. 2001; Siebenmorgen et al. 2004a) 
and theoretical (Voit 1992) studies. The latter suggest that PAH molecules hit by single 
energetic EUV/X-ray photons can be efficiently destroyed by photo-thermo dissociation or 
Coulomb explosion. Since AGN are copious emitters of hard photons, PAH molecules near 
AGN will be destroyed unless shielded by a large obscuring column. 

(2) The low excitation fine-structure emission lines like [Nell] 12.8/im are among the 
dominant emission lines of HII regions. Observations of starburst galaxies, as well as a 
combination of evolutionary synthesis and photoionization modelling, show [Nell] 12.8/im to 
be stronger than higher excitation mid-infrared lines ([NeV], [OIV], [Nelll]) in typical ionized 
regions excited by young stellar populations (Thornley et al. 2000; Verma et al. 2003). Use of 
low excitation lines as star formation tracers requires, however, the consideration of possible 
contributions from the AGN Narrow Line Region (NLR) which can be significant despite 
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the generally higher excitation of such regions (e.g. Spinogho & Malkan 1992; Alexander et 
al. 1999) compared to starburst H ll regions. 

Section 2 of the paper describes the sample, observations and data reduction used to 
obtain the line and continuum fluxes in our sources. Emission lines that are relevant to the 
present study are tabulated for all sources. In Section 3 we discuss the widespread presence 
of PAH emission and its relation to other components of the QSO spectra. Finally, Section 4 
addresses the issue of star formation in host galaxies of QSOs, shows the importance of this 
process and compares our results with earlier findings. In a forthcoming paper, we discuss 
the implications of our results for QSO SEDs in general. We adopt flm = 0.3, fl\ = 0.7 and 
Hq = 70 kms^i Mpc-^ 



2. The QUEST PG QSO sample: Observations and reduction 

2.1. The sample 

As part of the Spitzer spectroscopy project QUEST (PID 3187, PI Veilleux) we are 
studying QSOs, ultraluminous infrared galaxies, and the possible evolutionary connection 
between the two using the infrared spectrograph IRS (Houck et al. 2004). The QSO sam- 
ple is largely drawn from that of Guyon (2002) and Guyon et al. (2006). It consists of 
Palomar-Grccn (PG) QSOs (Schmidt & Green 1983) and covers the full ranges of bolomet- 
ric luminosity (~ 10^^'^^^'^Lq based on the B band absolute magnitude and the SED of Elvis 
et al. 1994), radio loudness, and infrared excess (\og{uL^{60fim) /Lboi)^ 0.02-0.35) spanned 
by the local members of the PG QSO sample (see also Jester et al. 2005, for a recent view 
on selection effects in the PG sample). B2 2201-I-31A is not a PG QSO but is in the sample 
because its B magnitude actually satisfies the PG QSO completeness criterion of Schmidt 
& Green (1983). It is one of the five radio-loud systems in our sample. At the sample's 
maximum redshift of 0.325, important emission lines like [OIV] 25.89//m stay within the 
IRS spectral range for all objects. The QUEST sample used in this paper includes 23 of 32 
objects from the Guyon sample. We exclude here two recently observed QUEST objects that 
are not yet fully processed by the Spitzer pipeline. We add two Palomar-Green objects from 
the Guyon et al. sample previously observed by Spitzer (PG0050+124 = IZwl; Weedman 
et al. (2005) and PG0157+001 = Mrkl014; Armus et al. (2004)) and two PG QSOs from 
another project (PID 20241, PI Lutz). Table 1 lists names and redshifts of all 27 QSOs in 
our sample. In total, our sample covers a range from -21 to -26, with median -23.3. In 
the remainder of this paper, we will compare some aspects of the PG QSOs to ultraluminous 
infrared galaxies whose properties will be presented in more detail in an upcoming paper 
based on QUEST data (Veilleux et al., in preparation). 
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2.2. Data reduction and line and continuum measurements 

For the QUEST sample, spectra were taken both at 5-14yum in the low resolution (SL 
short-low) mode and at 10-37/im in the high resolution (SH short-high and LH long-high) 
modes of the IRS . Slit widths of 3.6" to 11.1" include much of the QSO hosts as well as the 
vicinity of the AGN. Our data reduction starts with the two-dimensional basic calibrated 
data (BCD) products provided by version 12 of the Spitzer pipeline reduction. We used 
our own IDL-based tools for removing outlying values for individual pixels and for sky 
subtraction and SMART (Higdon et al. 2004) for extraction of the final spectra. Small 
multiplicative corrections were applied to stitch together the individual orders of the low 
resolution and high resolution spectra, as well as additive corrections for residual offsets 
still found between the low resolution spectra and the SH and LH high resolution spectra 
after zodiacal light correction of the latter. Emission line fluxes were measured using flts 
of Gaussian lines superposed on a local continuum. PAH fluxes were mostly measured by 
simultaneously fltting Lorentzians to the main 6.2, 7.7, and S.d/im features superposed on 
a 5.3 to 9.6//m (rest frame) continuum approximated by a second order polynomial, and in 
a few cases by combined Lorentzian and continum flts over smaller ranges. We use below 
the flux for the brightest (7.7;um) feature as PAH strength indicator Fpah- In PG0050-I-124, 
where this feature was present but difficult to quantify in the presence of strong silicate 
emission, we have estimated the flux scaling from the flux of the well-detected 6.2/im PAH 
feature, using the 7.7/6.2 feature ratio measured in the starburst/AGN NGC6240 (Lutz et 
al. 2003; Armus et al. 2006). Flux upper limits {3a) were derived adopting typical widths 
for the lines (~600km/s) and broad features (~ 0.6//m for the 7.7//m feature). In one 
source (PG 1307-1-085) we could only analyse the high resolution spectra since low resolution 
data are proprietary to another program. This limits our sample to 26 objects for analysis 
related to PAH emission. QUEST data for ultraluminous infrared galaxies that are used for 
comparison with the QSOs were processed in the same way. 

Figure 1 shows individual example spectra of two PG QSOs illustrating the broad differ- 
ences between QSOs with strong and weak PAH emission. Table 1 also fists the continuum 

flux densities in several mid- and far-infrared bands. In the mid-infrared, the quality of the 
IRS spectra is superior to most pre-existing IRAS or ISO pfiotometry. For that reason and 
to characterise the continuum at wavelengths reasonably free of emission features, we have 
derived average observed flux densities over narrow bands (AA/A ~ 0.07) centered at several 
rest wavelengths. The shortest wavelength continuum point at 6/i.m is shortward of the main 
PAH complex; the 15/xm point is between the two sificate emission peaks at about 10 and 
18/im (Siebenmorgen et al. 2005; Hao et al. 2005; Sturm et al. 2005), but still partly affected 
by silicate emission if present. The longest wavelength point at 30//m is near the long end 
of the IRS spectra and already beyond the strongest part of the longer wavelength sificate 
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feature. We later use continuum luminosities at these rest wavelengths: 6/xm, 15//m, and 
30 /im, defined as i/L^. 

Far-infrared fiuxes have been taken from the literature, usually giving preference to ISO 
based results over IRAS based ones, because the smaller effective ISO beams reduced the 
susceptibility to cirrus contamination at 100//m. The far IR luminosity, L(FIR), is often 
obtained from Ffir which is defined as FpiR = 1.26(Sioo + 2.58 x Sgo) x lO^^^W/cm^, where 
Sioo and Sgo are flux densities in Jy (e.g. Sanders & Mirabel 1996). In order to use a consistent 
definition over the full z<0.33 redshift range of our sample and to reduce the sensitivity to 
galactic cirrus contamination that is most problematic at lOOyum, we adopt a far-infrared 
luminosity z/Lj, which is solely based on the flux at rest wavelength 60/im interpolated from 
the photometry in the literature. Even nearby QSOs are typically close to the far-infrared 
detection limits of the previous space missions. In six cases where only one of Sioo or Seo was 
detected we estimated the other flux as the lower of the measured limit and an extrapolation 
using the detected flux and the ratio Sgo/Sioo = 0.93 that is the median for the part of the 
sample detected in both bands. The flux at 60/im rest wavelength was interpolated linearly 
between the observed 60 and 100/im points. For seven PG QSOs in our sample that are 
undetected in both bands, upper limits on uL^, at 60/im rest wavelength were derived in an 
analogous way interpolating linearly between the limits at observed wavelengths of 60 and 
100/xm. Since QSO fluxes are close to the detection hmits of the original IRAS and ISO 
references and because of residual cirrus contamination, uncertainties of order 30-50% for 
the far-infrared fluxes may well occur. 

3. Results 

Table 2 lists the intensities of the three emission lines used in the present analysis, 
[Neil] 12.8/im, [Nev] 14.3/im, [0iv]25.9/xm and of the emission feature PAH7.7/xm. The 
table includes detections as well as upper limits. A more comprehensive list including more 
lines will be given in a forthcoming paper. 

3.1. PAH emission is common in QSOs 

As is clear from inspection of Table 2, some of the QSOs in our sample show only high 
excitation fines, likely originating in the NLRs of these objects. Others (11 of 26) show 
prominent PAH 7.7//m and other PAH features, indicative of signiflcant star formation. The 
existence of "composite" sources (see S02 for deflnition and references to ISO-based studies) 
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that show both AGN and starburst properties, has been known for several years. Such sources 
contain both powerful star forming regions, giving rise to the observed PAH emission through 
near-UV heating, and high luminosity NLRs that are excited by the central radiation source. 
Since the star forming activity is a property of the host galaxy, and the NLR line luminosities 
scale with the AGN continuum luminosity, the intensity of the starburst emission features 
relative to the continuum might be expected to be weaker in AGN whose luminosities are 
in the QSO regime. Because the S02 sample included only three objects with QSO-hke 
luminosities (i.e. L(FIR)> 5 x 10^^ Lq) and suffered from the limited S/N achievable with 
ISO, a quantitative verification of this idea was not possible up to now. 

Our new sample includes 27 high luminosity AGN and it is therefore important that 
the fraction of clear PAH emitters in the sample is large. This is strengthened further by 
the detection of PAH 7.7//m and 11.3//m peaks in the average spectrum of the 15 QSOs not 
showing individual PAR detections (Fig. 2, bottom). Since this average spectrum excludes 
sources with individual PAH detections, the PAH emission seen is unlikely to be due to a few 
PAH-strong sources only, in an otherwise largely PAH-free group. It also has the implication 
that typically the true PAH fluxes of the PAH nondetections cannot be far below our limits. 
This high incidence of PAH emission is our first major conclusion and is further discussed 
in the following sections. 



3.2. Trends with level of PAH emission 

Before embarking on an analysis of the correlations of different starburst and AGN 
tracers, we use the two average spectra of the PG QSOs that are individually detected in 
PAH emission (11 objects) or individually undetected in PAH emission (15 objects) (Fig. 2) 
to identify some of the salient trends in our data. The two spectra have been obtained by 
averaging the individual spectra after normalizing to the same total flux in the mid-infrared 
5-25//m rest wavelength region. We caution that signiflcant variations are present also within 
those two groups. 

(i) PAH7.7/im is 2.5 times stronger (relative to the total 5-25/xm mid-infrared flux) in 
the average of the objects with individual PAH detections. This is a natural consequence 
of the object grouping for the two average spectra. As noted, the relatively small difference 
argues for the presence of PAH emission in most of the objects that did not have individual 
detections. 

(ii) The far-infrared continuum emission is relatively stronger compared to the mid- 
infrared continuum in the objects with PAH detections. This is apparent in several ways. 
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Firstly, all PAH detections are also far-infrared (60//m rest-frame) detections while 7 of 

15 PAH nondetections arc also not detected in the far infrared. Secondly, the mean ratio 
z/iS,(60/im)/z/F^(6/im) is 1.72±0.54 for the 11 PAH detections, 1.00±0.15 for the 8 FIR- 
detcctcd PAH nondetections, < 0.88 for the 7 objects undetected in both far-infrared and 
PAH, and finally <0.94 for all 8+7 PAH nondetections. The values for the two upper 
limits assume that the far-infrared fluxes are less than or equal to the measured limits for 
the individual far-infrared nondetections. Thirdly, inspection of the average IRS spectra 
in Fig. 2 shows significant differences in the extrapolation to 60//m: The average spectrum 
of the PAH detections continues to rise beyond 25//m, while there is an indication for a 
downturn in the average of the PAH nondetections beyond this wavelength. Regardless of 
uncertainties in the extrapolation that are related to the assumed intrinsic spectral shape 
and to technical issues like IRS spectral response calibration and zodiacal light subtraction, 
the 60;um to 6/xm flux density ratio must be signiflcantly larger in the average spectrum of 
the sources with detected PAHs. We estimate that this ratio is larger for the PAH detections 
by a factor 2 to 3. 

We have also tentatively grouped the PAH nondetections into two average spectra for 
the 8 far-infrared-detected and 7 far-infrared-nondetected sources. While the statistics and 
S/N are poorer in these samples, there are indications that the sources undetected in both 
tracers are at the end of the physical trends described in (i) and (ii) - little PAH, weak 
far-infrared, and indication for more of a downturn at A > 25//m; that is, they are not just 
fainter overall. 

(iii) The lower excitation [Nell] line is stronger in the QSOs detected in PAHs. Nor- 
mahzed by the higher excitation hues [SIV] and [OIV], the [Nell] line is 1.9 and 1.7 times 
stronger in the QSOs detected in PAHs compared to those undetected in PAHs. 

(iv) The broad 10/im silicate emission peak is apparently weaker in the sources with 
detected PAHs. We argue this is largely an artefact of having a stronger starburst component 
with its PAH emission. This is illustrated by the dotted hne in the upper panel of Fig. 2, 
which shows the average spectrum after subtracting a PAH-dominated M82 spectrum (Sturm 
et al. 2000) roughly scaled to the PAH features, and the dashed hne, which shows the 
equivalent result using the average of 12 starburst-dominated QUEST ULIRGs (see §4.1) 
as a starburst template. As noted by, e.g., Sturm et al. (2000), a main difference between 
different starburst templates is in the level of increase of the A > 12/xm very small grain 
continuum, even for similar shorter wavelength PAH spectra. Our two choices illustrate the 
effect of such a template variation. The two PAH complexes around 7 and 12/xm and the 
far-infrared upturn serve to 'fill in' the minima between and around the two silicate emission 
peaks. This is analogous to the well-known difficulty of quantifying silicate absorption in 
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spectra with strong PAH emission. We note that subtracting a starburst spectrum from 
the PAH-strong average QSO spectrum, in addition to recovering the correct pronounced 
sihcate emission, also indicates a flatter or decreasing extrapolation to beyond 25//m, similar 
to what is seen in the PAH-weak average QSO spectrum. 

(v) Molecular hydrogen emission in the S(l), S(2) and S(3) rotational hues has a larger 
equivalent width in the PAH-strong average spectrum, by factors 3-5 comparing these de- 
tections with the detections/limits for the PAH-weak average spectrum. 

All these trends are consistent with a starburst component (containing strong PAH and 
far-infrared continuum, low excitation fine-structure line, and possibly molecular hydrogen 
emission) being superimposed in increasing proportion on a pure AGN spectrum (consisting 
of warm and hot dust continuum, silicate, and high excitation line emission) . 

4. Discussion 

4.1. Nature of the QSO far-infrcired emission and the steirburst-AGN 

connection 

We use our sample to compare several AGN- and starburst related quantities observed 
in the PG QSOs, in an attempt to identify the likely origin of the far-infrared emission. We 
compare the starburst tracers in these QSOs with those of pure starbursts from two samples: 
(i) A subset of starburst-dominated QUEST ULIRGs observed with IRS ('SB-ULIRGS')- 
In order to restrict ourselves to the ULIRGs with the highest starburst contribution to 
their infrared luminosity, we require these objects to have no [0iv]25.9//m detection, a 
peak ratio of the PAH7.7/xm feature to its local continuum of at least 1, and a 5-10//m 
spectrum in which visual inspection shows that the absorption features like the Q^m. ice 
feature are not dominant, although sometimes present. The last criterion is used to avoid 
ambiguities concerning the internal energy sources of the most heavily obscured ULIRGs 
(Spoon et al. 2004a,b). (ii) A small sample of six /5'0-observed local starbursts (M82, 
NGC253, NGC1808, IC342, NGC3256, NGC7552) for which the PAH emission is measured 
without significant aperture corrections relative to the far-infrared data from IRAS, ISO or 
the Kuiper Airborne Observatory KAO. The ratio of PAH and far-infrared emission is known 
to vary somewhat with physical conditions (see e.g. discussion in §4.2 of Lutz et al. 2003), 
decreasing with the average intensity of the radiation field (e.g., Dale et al. 2001). Because of 
these trends with ISM conditions, we add as another group of comparison objects (iii) twelve 
'FIR quiescent' galaxies (with particularly low Lfjji/Lb and low Sgo/Sioo) from the sample 
of normal galaxies of Lu et al. (2003), for which we convert their PAH fluxes measured in 
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ISOPHOT-S spectra to our measurement procedure and apply aperture corrections based on 
their PAH-dominated ISOCAM LW2 images. We add to this group NGC 891 that has been 
mapped with ISOPHOT-S by Mattila et ah (1999), arriving at a sample of 13 FIR-quiescent 
objects. The three comparison samples cover a range in luminosity as well as in the average 
intensity of their interstellar radiation fields that is reflected in their mid-to far-infrared 
SEDs. 

An important new result is the clear correlation between L(PAH7.7/im) and L(FIR) for 
the PAH-containing QSOs. Figs. 3 and 4 compare the luminosities and fluxes of PAH7.7/im 
and 60/um continuum. The QSOs with PAH detections and the starburst-dominated ULlRGs 
follow the same trend. Fig. 5 shows the equivalent continuous trend from QSOs to starburst- 
dominated ULIRGs for the comparison of [Nell] 12.8/xm and 60/xm continuum. 

These trends also imply almost identical mean ratios for the QSO and starburst-dominated 
ULIRG populations. For the 11 QSOs with both PAH and far-infrared emission detected we 
find < L(PAH)/L(FIR) 0.0110 ± 0.0021 while for the 12 starburst-dominated ULIRGs 
we get < L(PAH)/L(FIR) >= 0.0130 ± 0.0015. Similarly, comparing [Neli] 12.8//m and far- 
infrared emission we get for the 18 QSOs with both quantities detected < L([NeII]) /L(FIR) >= 
(5.15 ±0.73) X 10-^ and for the 12 starburst-dominated ULIRGs < L([NeII])/L(FIR) >= 
(4.22 ±0.57) X 10~^. In both cases, the ratio of the star formation indicator and far-infrared 
emission is the same in QSOs and in starburst-dominated ULIRGs - the observations are 
consistent with starbursts producing all of the QSO far-infrared emission if their specific 
properties are similar to those in starburst dominated ULIRGs. 

This argument depends on the adopted specific properties of the comparison starforming 
galaxies, because of the changes of the PAH to 60/im ratio with ISM properties mentioned 
above. For ULIRG starbursts we obtained < L(PAH)/L(FIR) >= 0.0130 ± 0.0015. For the 
six lower luminosity starbursts for which the ratio of PAH and far-infrared can be derived 
from ISO data without significant aperture corrections, we obtain 

< L(PAH)/L(FIR) >= 0.0399 ± 0.0087. Finally, for the sample of 13 FIR-quiescent normal 
galaxies we get < L(PAH)/L(FIR) >= 0.126 ±0.016 with individual objects like NGC 891 
reaching up to a value 0.2. 

Comparing these three ratios to the average of the QSOs, we find that almost all, at least 
~30%, and at least ~10% of the QSO far-infrared emission would be due to non-AGN sources 
(that is current and recent star formation), adopting the values for ULIRGs, starbursts and 
FIR-quiescent galaxies as templates. Irrespective of the adopted template, PAH-based values 
are lower limits since for QSOs there is the interesting possibility of additional weakening of 
PAH emission in star-forming regions that are also exposed to AGN radiation, by destruction 
of the PAH molecules. This may be particularly relevant if some of the far-infrared emission 
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is generated by star formation occurring in compact clusters or disks very close to the AGN 
(Davies et al. 2004a,b). This makes values towards the upper end of the range for the 
non-AGN contribution more likely. 

Using only PAH and far-infrared evidence, a formal solution is possible which maximizes 
the AGN fraction of the far-infrared emission by assigning the PAH emission to a very FIR- 
quiescent galaxy. The QSO PAH luminosities are on average larger than those of the FIR 
quiescent galaxies, but the most luminous quiescent objects are of a similar PAH luminosity 
as the typical QSOs (Fig. 3). Comparing these objects to QSOs, the AGN would then 
contribute roughly 90% of the QSO far-infrared flux. These objects yield biased comparisons, 
however, since they are specifically selected for FIR quietness. A fairer comparison would 
need to populate the PAH/FIR plane of Fig. 3 with a large and complete sample of non-AGN 
galaxies. While such a sample is currently not available, we qualitatively indicate in Fig. 3 
the expected locus of such non-AGN galaxies by connecting the locations of the averages 
of our three small comparison samples (dotted line). While the FIR quiescent objects are 
likely more PAH luminous than the average galaxy of the same FIR luminosity, and thus 
their point placed too high, the overall trend of PAH flux with FIR flux is certainly robust, 
and the slope of this relation for non-AGN galaxies less than 1, in agreement with the 
photometric relation between total infrared and ISOCAM LW2 emission studied by Chary 
Sz Elbaz (2001). Applying such a luminosity-dependent ratio of PAH to FIR luminosity to 
QSOs is consistent with our finding of star formation dominating the FIR emission in many 
of our QSOs. The alternate scenario, in which host galaxies are all PAH luminous but FIR- 
quiescent and the far-infrared emission of QSOs is AGN dominated, is both complex and 
inconsistent with some of the evidence presented in §3.2. Specifically, it cannot explain why 
both the far-infrared and the PAH emission rise together relative to the mid-IR continuum 
(i.e. AGN) luminosity (see also Figure 2). Moreover, this scenario does not agree with the 
multiwavelength evidence for active star formation mentioned in §4.3. 

Considering these trends and the fact that the far-infrared luminosities of the PG QSOs 
are typically in the > IO^^Lq regime, we conclude that for the average QSO in the sample the 
star formation contribution to the FIR emission is at least 30% (applying the starburst PAH 
to 60/xm conversion), and that star formation may well be dominant. Comparing the locus 
of QSOs in Fig. 3 with the trend for star-forming comparison objects suggests a tendency 
for the star formation contribution to be largest in the most FIR-luminous QSOs. 

Similar considerations can be made for the comparison of [Nell] 12.8/im and far-infrared 
emission in QSOs and ULIRGs (Fig.5). Here, the sources of uncertainty are the NLR contri- 
bution to [Nell] 12.8/xm for the QSOs and the possibility of different extinction in the QSO 
starbursts and the ULIRGs which show considerable mid-IR extinction (Genzel et al. 1998). 
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The similarity between QSOs and starburst-dominated ULIRGs in trends and ratios 
based on PAH7.7/im, [Nell] 12.8/xm, and far-infrared continuum no longer holds when we 
plot a clearly AGN- related quantity on one axis. Fig. 6 shows that the ratio of PAH and 6fim 
continuum differs between QSOs and starburst-dominated ULIRGs by more than an order 
of magnitude, and that the two classes separate clearly in the diagram. An additional strong 
non-starburst component is required in the QSOs which clearly is the AGN-heated warm 
dust continuum. We note, however, a clear correlation between the luminosities of AGN 
6//m continuum and starburst PAH among the QSOs. We will now argue on the basis of a 
more comprehensive comparison that this correlation is indirect, caused by a star burst- AGN 
connection in QSOs. 

In a flux-limited sample selected based on an emission component that is directly due to 
the AGN, as is the case for the PG sample (selected by pointhke appearance, B magnitude 
and blue U-B color), we expect correlations between the luminosities of the various AGN 
tracers to arise directly as a consequence of the sample selection. Correlations are not ex- 
pected when comparing AGN- and starburst-related quantities unless there is a real physical 
correlation. 

A summary of correlation coefficients among the continua, PAH and line fluxes measured 
in the QSOs is given in Table 3. Some of the relations arc also shown in Figs. 3 to 8. In 
Table 3 we list Spearman rank correlation coefficients for the detected objects and their 
significance both for himinosities and for the observed fiuxes. Given the selection of the 
PG sample over a relatively narrow fiux range, the correlations in fiux arc generally less 
tight, but agree with the findings based on luminosities. Our dataset includes upper limits, 
in particular for feature measurements and far-infrared fluxes, as well as trends between 
luminosity and distance caused by the flux limited selection. For these reasons, we also list 
in Table 3 partial correlation coefficients (explicitly excluding the effect of distance) that 
have been computed using the formalism of Akritas & Siebert (1996) which provides partial 
correlation analysis for censored data. The results from this analysis of the full censored 
dataset are in agreement with those from the detections only. 

A flrst and expected flnding is the signiflcant correlations between various luminosities 
tracing the AGN components, e.g. 6/im continuum and the high excitation emission lines. 
There is a noticeable spread in the ratio of mid-infrared continuum emission and emission in 
the mid-infrared high excitation lines [Oiv] 25.9;um and [Nev] 14.3;um, however, equivalent 
to a spread in equivalent width of these lines. This is not uncommon in other AGN narrow 
emission lines like [O ill] A5007 (e.g. Boroson & Green 1992; Baskin & Laor 2005) and has 
implications for the reliability of NLR lines as direct tracers of AGN bolometric luminosity. 
Netzer et al. (2006) discuss this large spread of [O ill] A5007 equivalent width which is likely 
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a general property of AGN and also a function of source luminosity (the 'Baldwin-effect') 
and emission line reddening. 

Clear correlations are also seen between quantities tracing starburst activity (PAH l.T/im 
being the cleanest) and others that in the QSOs must be almost fully dominated by the AGN 
(like 6//m continuum or high excitation lines). Such correlations are not caused by the PG 
sample selection by B-band flux and U-B color, and must indicate a true relation of more 
luminous QSOs on average being associated with more luminous star formation. In the 
presence of such a "starburst- AGN connection", evidence on the causal links connecting 
the correlated obscrvablcs has to be obtained from the quality of the correlations, and in 
particular from comparing the absolute values of the correlated quantities with templates as 
done above for far-infrared and PAH. 

While the size of the present sample is modest. Table 3 shows that the correlation 
between PAH7.7/xm and far-infrared luminosities is one of the tighter among the combina- 
tions we investigated. Consistent with this trend is the stronger correlation of [Nell] 12.8/xm 
with far-infrared than is the case for far-infrared with either [Oiv] 25.9/im or [Ne v] 14.3/im. 
[Neil] 12.8//m is a hue that is emitted both in the NLR of AGN (Sturm et al. 2002), with 
higher excitation lines usually being stronger, and as the strongest mid-infrared line in the 
spectra of most starbursts (Thornlcy ct al. 2000; Vcrma ct al. 2003). The good corre- 
lation of [Nell] 12.8/um and far-infrared thus indicates a strong starburst contribution to 
both [Nell] 12.8/im and far-infrared emission, reinforcing the conclusion reached from the 
starburst-likc ratios of these quantities. In a larger sample, these considerations could be 
expanded to a more rigorous test on the basis of the quality of the correlations. Even in 
the presence of a starburst- AGN connection, starburst tracers should usually correlate more 
tightly with other starburst tracers than with AGN tracers. Our sample is not big enough 
for robust conclusions of this type. Looking at the probability of exceeding the partial cor- 
relation coefficients for the full censored dataset in the null hypothesis of uncorrelated data 
(column 11 of Table 3), it is nevertheless reassuring but certainly tentative that the seven 
least significant of the 14 correlations discussed arc 60/im vs. Gfira, 60/im vs. [Oiv] 25.9;um, 
60/im vs. [Nev] 14.3/im, PAH vs. 6/xm, PAH vs. [O iv] 25.9/im, PAH vs. [Ne v] 14.3;um, 
6/xm vs. [Ne v] 14.3/xm. With the exception of the last, these less significant correlations are 
all of the type starburst tracer vs. AGN tracer, in the superposition scenario outlined in 
§3.2, and provided that the far-infrared is counted as a star formation tracer. 

We conclude that, while there arc also 'indirect' correlations caused by a global correla- 
tion of AGN and starburst luminosity in our PG sample, the relation between PAH, [Nell], 
and far-infrared is real and reflects the starburst component. The most important support 
for this interpretation is the starburst-like ratios of these three quantities. Starbursts con- 
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tribute at least ~30% and likely most of the far-infrared emission in the average QSO in our 
sample. An upper limit to the starburst contribution is imposed by the need for a realistic 
continuation to longer wavelengths of the AGN mid-IR continuum, which cannot fall off 
more steeply than the Rayleigh- Jeans like emission of optically thin dust of an appropriate 
temperature. The true slope is likely somewhat shallower due to variation in temperature 
and due to non-negligible optical depth in part of the mid-IR emitting region. We feel that 
the origin and interplay of silicate emission and continuum in the SED of the AGN is not 
yet well enough measured or modelled for an accurate AGN continuum extrapolation of this 
type. This is also due to the need for unambiguous decomposition of starburst, silicate, and 
AGN continuum. Nevertheless, we consider the minimum pure AGN far-infrared continuum 
required by the data to be consistent with our global conclusion from the PAH7.7//m and 
[Nell] 12.8/xm emission. 



4.2. Mid-infrared diagnostics and the starburst- AGN connection in QSOs 

Extending earher ground-based work (e.g., Roche et al. 1991), Genzel et al. (1998) and 
Laurent et al. (2000) presented the first empirical versions of the tools now used to separate 
AGN-powered from starburst-powered infrared galaxies. The basis of these tools is that 
the intensity of the PAH features in starburst-powered systems traces the starburst's (far- 
infrared) luminosity, while these features arc easily destroyed by the strong and hard AGN 
radiation. Our present study widens the scope of the Genzel et al. (1998) work by focusing 
on star formation signatures in high luminosity bona-fide AGNs. We have used our QUEST 
QSO sample to look for three starburst signatures, strong PAH features, strong [Nell] lines 
and strong far- infrared continuum. Although two of those ([Nell] and far-infrared) can also 
partly originate in AGN environments, we have argued in §4.1 that the three quantities scale 
with each other and are tracing significant star formation in most objects of our PG QSOs 
sample. The measured far-infrared luminosity z/L,y(60/im) ranges between 1.7 x 10^'^ and 
2.5x10^^^0, covering a wide range of starburst luminosity up to the ULIRG regime, and 
limits for the remaining QSOs are consistent with starburst emission in the same range of 
luminosities. 

An important result is the correlation of PAH (starburst) luminosity and AGN luminos- 
ity in our sample. This extends to higher luminosity a similar result obtained by S02 on the 
basis of ISO spectroscopy of mostly lower luminosity Scyfcrts, and of a wide range of optical 
and near-infrared studies suggesting elevated starburst activity in Seyferts (e.g. Heckman 
et al. 1997; Oliva et al. 1999; Gonzalez- Delgado et al. 2001; Imanishi 2003; Kauffmann et 
al. 2003). Such a connection between small scale AGN feeding and larger scale starburst 
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activity is plausible (e.g., Norman & Scoville 1988) and may play a role in establishing the 
black hole mass to bulge velocity dispersion relation in galaxies. Its details are far from 
trivial, however, and warrant observations with higher spatial resolution to elucidate the 
spatial structure of star formation in these QSOs (see for example Cresci et al. 2004). 

Another effect of this connection relates to the interpretation of correlations between 
QSO properties mcasiu^cd at different wavelengths. Some observed correlations may be 
indirect, driven by the starburst-AGN connection, as argued for example by Rowan- Robinson 
(1995) for optical and far- infrared continua. Haas et al. (2003), however, have used among 
other arguments the observed correlation between rest frame mid-infrared and far-infrared 
continuum in their QSO sample to argue for an AGN origin of the latter. The sensitivity of 
the ISO spectroscopic data they had available did not allow for a conclusive test on the basis 
of PAH emission. While we confirm the mid- to far-infrared correlation for QSOs for our 
PG sample (e.g.. Fig. 7), we use our higher sensitivity 5pzfeer PAH 7.7/im and [Nell] 12.8//m 
data to argue for an indirect nature of this correlation for the luminosity range covered by 
our sample, induced by a starburst-AGN connection. A similar test remains to be done for 
the highest luminosity members of the Haas et al. (2003) sample, which are not sufficiently 
represented in our local PG sample or in the IRS spectra of radio galaxies and radio-loud 
QSOs of Haas et al. (2005). Such high quality mid- infrared spectra of highest luminosity 
QSOs should also be able to test whether a trend for decreasing FIR to MIR ratio at highest 
optical luminosity (as suggested by Fig. 4 of Haas et al. 2003) reflects an increase in relative 
AGN intensity compared to the host and its star formation. 

A central question in the study of QSOs and ULIRGs is their possible evolutionary 
relation (e.g. Sanders et al. 1988). Our finding of luminous starburst activity in many QSOs is 
clearly consistent with such an evolutionary link between ULIRGs (which show ultraluminous 
star formation and frequent coexisting AGN) and QSOs (which show ultraluminous AGN 
and frequent coexisting starbursts). From such basically energy- related considerations, an 
evolutionary path with a clear time arrow is, however, difficult to demonstrate and distinguish 
from more random processes. Including structural and dynamical information, e.g. from 
other elements of the QUEST program, will better probe this link. 

4.3. Direct AGN heating of cold dust and PAHs? 

Several models have proposed a direct AGN heating of the far-infrared emission of QSOs 
(e.g., Sanders et al. 1989). A basic feature of such models, needed in order to fit QSO SEDs 
with moderately strong far-infrared emission, is a significant covering factor by obscuring 
dust at relatively large (few kpc) distances from the central AGN that is not shadowed by 
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matter closer in. In the model of Sanders et al. (1989), for example, this is accomplished 
by invoking a dusty galactic disk that is warped into the unshielded AGN radiation on such 
scales. In such a scenario, the star formation activity and associated PAH emission would be 
low. The PAH to far-infrared correlation that we find would have to be due to PAH excitation 
by the AGN itself at this relatively large distance. The required emission is significant - as 
we argued above, the ratio of PAH and far-infrared is similar to that in starburst-dominated 
systems, and far-infrared is a significant fraction of the bolometric luminosity for our sample 
(~10% mean for the systems with far- infrared detections). 

A main argument against this scenario is the likely destruction mechanism of PAH 
by AGN radiation: If this process works through destruction by individual EUV and X-ray 
photons, then the PAH carriers cannot survive at even kpc distances unless shielded by a large 
column able to stop the deeply penetrating hard photons (Voit 1992). Such a large absorbing 
column (Nh > lO^^cm"^) would, however, prevent the heating of a significant far-infrared 
emitting dust component at large distance, by absorbing the UV 'big blue bump' bulk part 
of the AGN SED. It would also absorb the near-UV (<13.6eV) radiation needed to actually 
excite infrared feature emission from the shielded PAH molecules. As argued by Voit (1992) 
and Maloney (1999), PAH molecules can survive even relatively close to a powerful AGN if 
placed behind large obscuration. Such high obscuring columns could be plausibly identified 
with the anisotropic obscuring structure postulated by unified AGN models. Exciting these 
surviving PAH molecules will however require a separate near-UV source, for example by 
reintroducing a circumnuclear star burst. An indirect transport of AGN near-UV radiation 
to this shielded material, for example by scattering of UV emerging into the AGN ionisation 
cones, appears unlikely to produce the required large PAH luminosities for ~1% scattering 
efficiency, as often assumed on the basis of polarimetric AGN studies (e.g. Pier et al. 1994). 

Preudhng et al. (2003) and Siebenmorgen et al. (2004b) present observations and models 
of radio-loud AGN in which they ascribe the full infrared SED including sometimes present 
PAH emission to the AGN, without invoking additional star formation. Their models (see 
discussion in Siebenmorgen et al. 2004a) invoke a difi^crent treatment of PAH destruction 
in the AGN radiation field. While this issue certainly deserves future study, their models 
may underestimate PAH destruction by the AGN. They predict strong PAH emission from 
optically thin dust illuminated by an AGN at a typical NLR distance, as well as significant 
PAH emission from an optically thick model for the central region of a nearby AGN which 
has most of its mid-infrared emission arising in a compact (~10 pc) region (Fig. 20, 21, 22 
of Siebenmorgen et al. 2004a). In contrast, spatially resolved observations put very strong 
limits on the PAH emission from such regions of some nearby AGN (Le Floc'h et al. 2001; 
Siebenmorgen et al. 2004a; Weedman et al. 2005; Mason et al. 2006). 
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Our interpretation is also consistent with other, partly circumstantial, evidence for star 
formation in QSOs. Canalize & Stockton (2001) find optical spectroscopic evidence for 
relatively recent < 300Myr star formation in a sample that ranges from ULlRGs to some 
of the moderately FIR-bright PG QSOs of our sample. Veilleux et al. (2006) find some 
PG QSOs brighter than the H-band fundamental plane, possibly indicating circumnuclear 
star formation. Molecular gas detections of PG QSOs (e.g. Evans et al. 2001; Scoville et 
al. 2003), while not directly probing star formation, suggest sufficient material to power the 
far- infrared emission if star formation is efficient. 



4.4. Compcirison to QSO steir formation estimates based on [O ii] A3727 

Our finding that star formation activity is able to power the far-infrared emission is 
in contrast to the result of Ho (2005). For a PG-based sample partly overlapping with our 
sample, he found insufficient star formation from an analysis using the [O ll] A3727 line as a 
star formation tracer. Large and uncertain corrections for extinction of this tracer are likely 
the main contributor to this discrepancy. 

Ho (2005) used [O ll] A3727 fluxes for his QSOs from the literature and the calibration of 
Kewley et al. (2004) for extinction- corrected [O ll] A3727 as a star formation indicator that is 
also a function of metallicity: SFR([OII],Z). Kewley et al. (2004) base their work on a sample 
of nearby field galaxies with detailed integrated optical spectroscopy, and star formation 
rates spanning 4 orders of magnitude centered on ~ IMQyr"^. They demonstrate that 
[O ll] A3727 is a star formation indicator as good as Ha over this range, provided extinction 
can be corrected for individual objects, and individual metallicities are known. 

To apply this cahbration to QSOs whose optical spectra are strongly dominated by 

AGN emission. Ho (2005) made three assumptions: (1) Screen attenuation of Av=l towards 
the star forming regions, (2) metallicity twice solar, (3) one third of the [O II] A3727 emis- 
sion comes from star formation. Another implicit assumption is (4) aperture corrections 
to [Oll]A3727 can be ignored for the hosts of this QSO sample with median z~0.09 ob- 
served with 2-5" optical spectroscopic apertures. Under these assumptions. Ho (2005) infers 
star formation rates of at most 20 MqJt~^ but often much less, and typically an order of 
magnitude below those inferred by ascribing the far-infrared emission to star formation. 

While all four assumptions contribute to the uncertainties of this approach, we believe 
that assumption (1) about the reddening to the star forming regions is the key to the sys- 
tematic difference to our Spitzer results. As noted, e.g., by Kewley et al. (2004), extinction is 
a strong function of intrinsic [O ll] A3727 luminosity. At the short wavelength of [O ll] A3727 
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where ~ 1.5 x Ay, increased dust extinction can relatively easily offset much of an in- 
crease in intrinsic line luminosity. Since the far-infrared luminosities of our quasars reach 
beyond lO^^L©, it is instructive to compare to the systematic optical spectroscopic study 
of a complete sample of luminous infrared galaxies by Poggianti & Wu (2000), which is ad- 
dressing objects with Lpjji ~ lO^^'^L© for our cosmology. More than half of their objects 
have so-called e(a) optical spectra, characterized by absorption in the higher Balmer lines, 
weak [On] A3727, and significant obscuration of the emission hues. E(B-V) is ~1.11 even in 
the simplifying screen assumption, which still ignores regions in such highly dusty objects 
that are too obscured to contribute to the optical hues. With observed L(Hq;)~ lO^^erg 
s~^ and [On] A3727/Hq; ~0.23, they are in the regime of or below the portion (1/3) of the 
QSO [O ll] A3727 luminosities that Ho (2005) ascribes to star formation, already considering 
that [O II] A3727 aperture corrections may be more significant for the Poggianti & Wu (2000) 
sample which has median z=0.0324 and similar spectroscopic apertures. Similar arguments 
apply, to a lesser degree, to the Poggianti & Wu (2000) objects with spectral classifications 
other than e(a). There is a large population of dusty luminous starbursts that could fit the 
weak [O ii] A3727 emission of the Ho (2005) QSOs. 

In summary, while the QSO optical data are consistent with the low star- format ion rate 
interpretation of Ho (2005), they are also consistent with much higher star formation rates. 
This is because the optical analysis of Ho (2005) is based on a single hue that is strongly 
extinction sensitive, and no extinction constraints are available for this component from 
the optical data. Less extinction sensitive data like the Spitzer infrared spectra are needed 
to break this degeneracy, and suggest more substantial star formation rates. Some of the 
QSO host optical spectra may share properties of infrared galaxies with e(a) optical spectra: 
little obscuration towards the stellar continuum spectra dominated by an older post-burst 
component, but still significant obscuration of the active star forming regions (Poggianti & 
Wu 2000). 



4.5. Implications for high redshift QSOs 

Deep submm and mm photometry has led to the detection of rest frame submm and 
far-infrared dust emission from radio-quiet QSOs at redshifts up to 6.42 (e.g. Omont ct al. 
2001; Isaak et al. 2002; Bertoldi et al. 2003). It has been variously argued through indirect 
arguments like CO content that this emission is star formation powered (e.g. Walter et al. 
2003), implying that these QSOs coexist with extremely powerful > lO^^L© starbursts. Our 
results support a starburst origin of QSO far-infrared emission, but do not extend to this 
luminosity and redshift regime. If the ratio of PAH to far-infrared emission for these QSOs 
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is similar to the one in the local QSOs and in ULIRGs, detection of PAH emission on top of 

a strong continuum may be within reach of Spitzer spectroscopy. PAH emission from similar 
luminosity SMGs is detectable at ULlRG-like ratios to the far-infrared emission (Lutz et 
al. 2005). Luminous PAH emission has also been reported in mid-infrared selected samples 
of z~2 infrared galaxes (Yan et al. 2005). Probing for PAH emission may currently be the 
only way to verify the assumption of simultaneous strong star formation and QSO activity 
in high redshift QSOs. 



5. Conclusions 

Sensitive Spitzer mid-infrared spectroscopy reveals the widespread presence of aromatic 
'PAH' emission features in z< 0.3 QSOs from the Palomar-Grecn sample, indicating the 
presence of powerful (i/L,^(60/im)~ 1.7 x 10^° to 2.5xlO^^L0) star formation activity in these 
systems. Starburst and AGN activity are connected in QSOs up to these high luminosities. 
By comparing the ratios of PAH7.7/im, [Nell] 12.8//m, and far-infrared emission in QSOs 
with starbursts we conclude that for the average QSO in our sample at least 30% and likely 
most of the QSO far-infrared emission is due to star formation. The data suggest a trend 
with the star formation contribution being the largest in the most FIR-luminous QSOs. 
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Note. — Col. (1) — Source name. 
Col. (2) — Rcdshift. 

Cols. (3-5) — Observed narrow band continuum flux densities around 6, 15, and 30 /xm rest wavelength, extracted 
from the IRS spectra. Missing entries are outside the rest wavelength range available for a given object. 
Cols. (6-8) — ISO or ffi^S' fluxes at observed 60 and 100 fim, and related references: FSC IRAS Faint Source 
Catalog, HOO - Haas et al. (2000), H03 - Haas et al. (2003), N86 - Neugebauer et al. (1986), S89 - Sanders et al. 
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(1989), SP - from IRAS Addscan/Scanpi 

Col. (9) — Luminosity distance in Mpc for a Ho=70kms~^ Mpc~^, fi^ = 0.3 and ^^A = 0.7 cosmology 
Col. (10) — Far infrared luminosity. See text for treatment of limits. 
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Table 2. Measured emission features 
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Table 3. Correlation of measured properties 
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Note. — Col. (1) — First variable. Wavelengths stand for continuum at that rest wavelength. 
Col. (2) — Second variable. 

Col. (3) — Number of sources detected in both quantities. 

Col. (4) — Spearman's rank correlation coefficient for luminosities, for detected sources. 



Col. (5) — Probability of exceeding the measured correlation coefficient for luminosities in the null hypothesis of uncorrelated 
data. Smaller values indicate more significant correlations. 
Col. (6) — Spearman's rank correlation coefficient for fiuxes. 

Col. (7) — Probability of exceeding correlation coefficient for fluxes in the null hypothesis. 

Col. (8) — Dispersion of Logio(A/B) for the detected sources. 

Col. (9) — Total number of sources measured (detections and limits). 

Col. (10) — Partial Kendall r coefficient describing correlation between luminosities A and B, excluding the effect of distance. 
Value computed using the formalism of Akritas &: Siebert (1996) that extends the Kendall r-coefficient to partial correlation in 
the presence of censored data (i.e. data including upper limits). 

Col. (11) — Probability of exceeding partial correlation coefficient in the null hypothesis. 
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Fig. 1. — IRS spectra of two PG QSOs, combined with far-infrared photometric fluxes. 
These two nearby sources were selected to illustrate the full range between sources with 
strong PAH, [Nell], and far- infrared emission and others with weak emission in all these 
tracers. The spectrum of PG1440-I-356 has been multiplied by a factor 4 for clarity. 
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Fig. 2. — Average spectra of QSOs. Individual spectra have been normalized to the same 
total 5-25//m flux before averaging. Note the change in spectral resolution from low at short 
wavelengths to high at > 9fim rest wavelength. Top: 11 QSOs for which PAH7.7;um is 
detected in the individual spectra. The dotted line is the same spectrum after subtract- 
ing an M82 spectrum scaled to the PAH features, and the dashed hne after subtracting a 
starburst-dominated ULIRG spectrum scaled to the PAH features. Emission hne residuals 
have been removed for clarity. The inset repeats the same spectra with different plot scaling, 
emphasizing the wider range SED trends. Bottom: 15 QSOs for which PAH 7.7//m is not 
detected individually. A broad maximum near 7.7/um rest wavelength as well as an 11.3/xm 
feature is detected in the average, however, indicating a high incidence of PAH emission in 
the contributing spectra. 
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Fig. 3. — Luminosities L(PAH 7.7/im) vs. L(60/im) for the QSOs. Starburst dominated 
ULIRGs, lower luminosity starbursts, and FIR-quiescent galaxiess are added for comparison. 
The thin dotted line connects the mean locations of these three groups of comparison objects. 
The small blue symbols repeat the QSOs without individual PAH detections, but assuming 
that the individual ratios of PAH7.7/im to rest frame 5-25/im flux are the same as in the 
average spectrum of Fig. 2 (bottom). A 30% la uncertainty is indicated in the lower right 
corner. 
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Fig. 4. — Fluxes F(PAH 7.7/im) vs. F(60/im) for the QSOs and for starburst dominated 
ULIRGs. 



- 36 - 



10 



10 



10' 



S 10° 

CN 



10^ \r 



10' 



1 1 llll 

- 
- 

^ QSOs 


Ill 1 1 1 1 1 1 1 1 




' 1 ' : 
- 
- 


_ + SB ULIRGs 






- 




+ 










+ 






^ ++_ 


+ 






















)K <-+ 














































































1 


III 1 





,1 



10 



10 



,1 1 



10 10 
uL, (60/^nn) [Lq] 



12 



10 



13 



Fig. 5. — Luminosities L([Nell] 12.8/im) vs. L(60yum) for the QSOs and for starburst domi- 
nated ULIRGs. 
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Fig. 6. — Luminosities L(PAH 7.7/im) vs. L(6/im) for the QSOs and for starburst dominated 
ULIRGs. 
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Fig. 7. — Luminosities L(6/im) vs. L(60/im) for the QSOs and for starburst dominated 
ULIRGs. 
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Fig. 8. — Luminosities L([0iv] 25.9/xm) vs. L(60/im) for the QSOs 



